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Abstract. We present the results from the dynamical analysis of the cluster of galaxies RXJ0152. 7-1357, which shows a 
complex structure in its X-ray emission, with two major clumps in the central region, and a third clump in the Eastern region. 
Our analysis is based on redshift data for 187 galaxies. We find that RXJ0152. 7-1357 appears as a well isolated peak in the 
redshift space at z = 0.836, which includes 95 galaxies recognized as cluster members. We compute the line-of-sight velocity 
dispersion of galaxies, cr v = 1322+gjj km s -1 , which is significantly larger than what is expected in the case of a relaxed cluster 
with an observed X-ray temperature of 5-6 keV. We find evidence that this cluster is far from dynamical equilibrium, as shown 
by the non Gaussianity of the velocity distribution, the presence of a velocity gradient and significant substructure. Our analysis 
shows that the high value of cr v is due to the complex structure of RX JO 152.7- 1357, i.e. to the presence of three galaxy clumps 
of different mean velocity. Using optical data we detect a low-velocity clump (with o" v = 300-500 km s -1 ) in the central South- 
West region and a high-velocity clump (with <x v ~ 700 km s -1 ) in the Eastern region, well corresponding to the South- West 
and East peaks detected in the X-ray emission. The central North-East X-ray peak is associated to the main galaxy structure 
with a velocity which is intermediate between those of the other two clumps and cr v ~ 900 km s -1 . The mass of the whole 
system within 2 Mpc is estimated to lie in the range (1.2-2.2)xl0 15 M G , depending on the model adopted to describe the cluster 
dynamics. Such values are comparable to those of very massive clusters at lower redshifts. Analytic calculations based on the 
two-body model indicate that the system is most likely bound, currently undergoing merging. In particular, we suggest that the 
South- West clump is not a small group, but rather the dense cluster-core of a massive cluster, likely destined to survive tidal 
disruption during the merger. 

Key words. Galaxies: clusters: general - Galaxies: clusters: individual: RX J0152. 7-1357 - Galaxies: distances and redshifts 
- Cosmology: observations 



1. Introduction lution is shown by the cluster population out to z ~ 0.3-0.4 

J-j ' , „ , . .., , „ , , „ (Adami et al. EOOOl Girardi & Mezzetti I200T1 but cf. Plionis 

rrt Clusters or galaxies are visible tracers or the network or matter u n mi „ , , rvci i_ 

. . i ° 20021. On the other hand, z > 0.5 clusters have more X-ray 

in the Universe, marking the high-density regions where fila- , . , . , T ,, . , i ^,-.,-.,! , 

° ° , ° substructures than lower-z clusters (Jeltema et al. 2005) and 

ments of dark matter loin together. In the hierarchical scenario . , . ,. c , . no , , . , , 

J b most clusters identified at z > 0.8 show an elongated, clumpy, 
of large-scale structure, clusters form via merging of smaller , cl . . ~ , ^ , . , | lrir , | 

b b b or possibly filamentary structure (e.g., Donahue et al. 1998 

clumps and accretion of material from large scale filaments „. . . . rnHin ..issmi., " . . . , 

b . Gioiaetal. 1999 Rosati 2004 1 thus suggesting that present ob- 
(e.g., Borgani & Guzzo 2001; Evrard & Gioia 2002 1. From , . , „ , . „ „ 

, „ servations are approaching the epoch of cluster formation. Our 

the observational side, signatures of past merging processes are u nvm ,,„_ . « o jj * ^ • 

b r t> bf results on RXJ0152. 7-1357 atz ~ 0.8 add further insights on 

found in cluster substructure and evidences for ongoing cluster ... . 

b b this issue, 
mergers are rapidly accumulating (e.g., Bohringer & Schuecker 

|2TK}2l BuoteHnOl Girardi & BivianoHnOl Evrard 120031 . The galaxy cluster RX J0152.7-1357 was disc overed in the 

Over the last few years significant progress has been made ROSAT Deep Cluster Survey (RDCS, Rosati et all 199811 in the 

to extend the above studies from local to distant clusters. ROSAT PSPC field rp60000rn00 observed in January 1992. 

Pioneering analyses suggest that no evidence of dynamical evo- 11 was independently discovered in the Wide Angle ROSAT 

Pointed Survey (WARPS, Ebeling et al. l2000> and reported in 

Send offprint requests to: M. Girardi (girardi@ts.astro.it) the Bright SHARC survey (Romer et al. 120001 . It appeared also 

* Based in part on observations carried out at the European Southern in the list of x _ ray extended sources obtained from Einstein 

Observatory using the ESO Very Large Telescope on Cerro Paranal Tpc data . 0ppenheimer et al flggfl . 

(ESO programs 166.A-0701, 69.A-0683, and 72.A-0759) and the ESO 

New Technology Telescope on Cerro La Silla (ESO program 61. A- The BeppoSax observations were used to derive a clus 
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erg s 1 (h — 0.5 and qo = 0.5), and a gas temperature kT = 
6.46+1$ keV (Delia Ceca et al. 12000b . RX J0152.7-1357 is 
characterized by a complex morphology with at least two 
cores, both in the optical and X-ray data as recovered by 
Keck imaging and Beppo-SAX data (Delia Ceca et al. 2000 1. 
Observations with Chandra also show a complex structure 
in the intra-cluster medium with the presence in the central 
cluster region of two peaks in the X-ray emission 95" apart 
(North-East: R.A.=l h 52 m 44 ? 18, Dec.=-13°57'15 ?84; South- 
West: R.A.=l h 52 m 39!89, Dec.=-13°58'27 ?48 [J2000.0]), 
and a possible third peak to the East (R.A.=l h 52 m 52 ?42, 
Dec.=-13°58'5:'52 [J2000.0]), see Maughan et al. ( 1277031 . 
The existence of an Eastern peak was confirmed by spec- 
troscopic VLT data and an independent analysis of the 
Chandra data by Demarco et al. 1 2005 1, who detect it at the 
> 3cr c.l. in X-rays (see their Fig. 1). Chandra observations 
gave a gas temperature for the North-East and South- West cen- 
tral X-ray clumps of kT = 5 .5 + °;^ k eV and kT = 5.2 + ' ' keV, 
respectively (Maughan et al. 2003 1. A complex structure with 
several clumps is also shown by the gravitational lensing anal- 
ysis of Jee et al. ©05 1: in particular, the mass clump A cor- 
responds to the Eastern X-ray peak. 

A number of evidences suggest that RXJ0152. 7-1357 may 
be undergoing a merger: the displacement between peaks of gas 
distribution and of galaxy/dark matter distribution (Maughan et 
al. l2003l Jee et al. 2005 1; the possible presence of a shock front 
(Maughan et al. 2003 ); the presence of galaxies showing a very 
recent star formation episode (J0rgensen et al. 2005 1; the segre- 
gation of star-forming and non star-forming galaxies probably 
induced by the intra-cluster medium interaction (Homeier et 
al. l2005ll . 

Demarco et al. (2005 1 have performed an extensive spec- 
troscopic survey of RX J0152.7-1357based on observations 
carried out with FORS1 and FORS2 on the ESO Very Large 
Telescope, obtaining more than 200 redshifts in the cluster 
field. Their analysis shows that RX J0152. 7-1357 is character- 
ized by a large velocity dispersion, ~ 1600 km s -1 , and indi- 
cates a very complex structure. In particular, the galaxy popula- 
tions inhabiting the regions around the three main X-ray peaks 
are characterized by different kinematical behaviour, in agree- 
ment with a cluster merging scenario. 

On the basis of Demarco et al. data we further investigate 
the internal dynamics of RXJ0152. 7-1357. The spatial and 
kinematical analysis of member galaxies is a powerful way 
to detect and measure the amount of substructure, to identify 
and analyze possible pre-merging clumps or merger remnants 
(Girardi & Biviano 2002 and refs. therein). This optical infor- 
mation is complementary to X-ray information since galax- 
ies and intra-cluster gas react on different time scales during 
a merger (see, e.g., numerical simulations by Roettiger et al. 
[19971 Ricker & Sarazin l200Tl Schindler l2002l 

The paper is organized as follows. We describe mem- 
ber selection and present our results for global proper- 
ties of RX J0152.7-1357 in Sect. 2. We present our anal- 
ysis of internal dynamics in Sect. 3. We discuss our re- 
sults suggesting a tentative picture of the dynamical status 
of RX J0152. 7-1357 in Sect. 4. We summarize our results in 
Sect. 5. 



Unless otherwise stated, we give errors at the 68% confi- 
dence level (hereafter c.l.) 

Throughout the paper, we assume a flat cosmology with 
Q. m = 0.3, Q A = 0.7 and H = 70 km s" 1 Mpc" 1 . For this 
cosmological model 1 arcmin corresponds to 458 kpc at the 
cluster redshift. 

2. Member selection and global properties 

Our data sample consists of the spectroscopic survey of 
RXJ0152.7-1357 presented by Demarco et al. J2005t . i.e. 187 
galaxies with available redshift (see their Tables 4 and 5). We 
assume a typical redshift error of 8 x 10~ 4 according to the 
authors prescriptions. 

The identification of cluster members proceeds in two 
steps, following a procedure already used for nearby and 
medium-redshift clusters (Fadda et al. 1996; Girardi et al. 
1996) Girardi & Mezzetti l20lHV 

First, we perform the cluster-member selection in velocity 
space by using only redshift information. We apply the adap- 
tive kernel method (Pisani 1993 ) to find the significant (> 99% 
c.l.) peaks in the velocity distribution. This procedure detects 
RXJ0152. 7-1357 as a well isolated peak at z = 0.836 assign- 
ing 103 galaxies considered as candidate cluster members (see 
Fig-D- Out of non-member galaxies, 61 and 23 are foreground 
and background galaxies, respectively. In particular, a second 
significant peak of 31 galaxies is shown at z — 0.638 suggest- 
ing the presence of a foreground system. 

All the galaxies assigned to the RX J0152. 7-1357 peak are 
analyzed in the second step, which uses the combination of po- 
sition and velocity information. We apply the procedure of the 
"shifting gapper" by Fadda et al. (1996 1. This procedure re- 
jects galaxies that are too far in velocity from the main body 
of galaxies and within a fixed bin that shifts along the distance 
from the cluster center. The procedure is iterated until the num- 
ber of cluster members converges to a stable value. We use a 
gap of 1000 km s _1 - in the cluster rest-frame - and a bin of 
0.6 Mpc, or large enough to include 15 galaxies. As for the cen- 
ter we consider the position of the biweight center, i.e. we per- 
form the biweight mean-estimator (ROSTAT package; Beers et 
al. [1990 1 for ascension and declination separately: this center 
is positioned between the North-East and South- West X-ray 
peaks (see § 1). The choice of using either one of the two X- 
ray peaks as cluster center does not affect the final results. 

The shifting-gapper procedure rejects eight galaxies to give 
95 fiducial members. The list of selected members corresponds 
to that in Table 4 of Demarco et al. ( 2005 1, but excluding galax- 
ies #306,509,557,650,895,1146,1239. Fig.|2]shows the plot of 
rest-frame velocity Vrf = (cz - (cz))/(l + (z)) vs. clustercentric 
distance R of galaxies in the main redshift peak. Finally, we 
recompute the biweight center on the 95 cluster members ob- 
taining: R.A.=l h 52 m 41 ?669, Dec.=-13° 57' 58:32 (J2000.0). 
Unless otherwise stated, we adopt this as cluster center. 

By applying the biweight estimator to cluster members 
(Beers et al. 1990 1, we compute a mean cluster redshift of 
(z) = 0.8357 + 0.0005. We estimate the line-of-sight (LOS) 
velocity dispersion, <x v , by using the biweight estimator and 
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Fig. 1. The redshift galaxy density, as provided by the 
adaptive-kernel reconstruction method. Unit on the y-axis is 
normalized to the density of the highest peak 

applying the cosmological correction and the standard cor- 
rection for velocity errors (Danese et al. 11980b . We obtain 
cr v = 1322^ km s , where errors are estimated through a 
bootstrap technique. 

To evaluate the robustness of the cr v estimate we analyze 
the integral velocity dispersion profile (Fig. [3}- The value of 
0" V (< R) sharply varies in the internal cluster region. A similar 
behaviour is shown by the mean velocity (V(< R)) suggesting 
that a mix of clumps at different redshifts is the likely cause 
for the high value of the velocity dispersion rather than indi- 
vidual contaminating field-galaxies. A robust value of cr v is 
reached in the external cluster regions where the profile flat- 
tens, as found for most nearby clusters (e.g., Fadda et al. 1996 1. 

The question of the presence of substructure is deferred to 
the following sections. Here we assume that the system is in 
dynamical equilibrium to compute virial global quantities. 

Following the prescriptions of Girardi & Mezzetti (2001 1, 
we assume for the radius of the quasi-virialized region R v ; r = 

0. 17 x o-y/H{z) = 2.0 Mpc (see their eq. 1 after introducing 
the scaling with H(z), see also eq. 8 of Carlberg et al. U997l for 
^20o)- Thus the cluster is sampled out to a significant region, 

1. e. R oat = 0.82 x R Vlr . 



We compute the virial mass (Limber & Mathews 1960J see 
also, e.g., Girardi et al. 1998 1 using the data for the N g observed 
galaxies: 



M = 3n/2 ■ crtRpy/G - C, 



(1) 



where C is the surface term correction (The & White 1986 1, 
and Rpy, equal to two times the (projected) harmonic radius, 



Fig. 2. Galaxies in the main peak of Fig.^ Left panel: rest- 
frame velocity vs. projected clustercentric distance; the appli- 
cation of the "shifting gapper" method rejects the galaxies in- 
dicated by open squares. Right panel: velocity distribution of 
all and member galaxies (dotted and solid histograms, respec- 
tively). 
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Fig. 3. Integrated mean velocity and LOS velocity-dispersion 
profiles (upper and lower panel, respectively), where (V) and 
cry at a given (projected) radius from the cluster center is esti- 
mated by considering all galaxies within that radius. The error 
bands at the 68% c.l. are shown. In the lower panel, the horizon- 
tal line represent X-ray temperature with the respective errors 
transformed in cr v imposing y6 spe c = 1 (see Sect. 4). 



Rpv = N g (N g - DK^jR;/), 



(2) 
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where Rjj is the projected distance between two galaxies. 

The estimate of cry is generally robust when com- 
puted within a large cluster region (see Fig. [3] for 
RXJ0152.7-1357 and Fadda et al. 173951 for other examples). 
The value of Rpy depends on the size of the sampled re- 
gion and possibly on the quality of the spatial sampling (e.g., 
whether the cluster is uniformly sampled or not). Here we ob- 
tain Rpy = (1.45 + 0.05) Mpc, where the error is obtained 
via a jacknife procedure. The value of C strongly depends 
on the radial component of the velocity dispersion at the ra- 
dius of the sampled region and could be obtained by ana- 
lyzing the velocity-dispersion profile, although this procedure 
would require several hundreds of galaxies. We apply the cor- 
rection obtained in the literature by combining data of many 
clusters sampled out to about R v [ r (C/My ~ 20%, Carlberg 
et al. [T9971 Girardi et al. [T9981 . We obtain M(< R out = 
1.65 Mpc) = (2.2 + 0.3) x 1O 15 M . Calling into question the 
quality of the spatial sampling, one could use an alternative 
estimate of Rpy on the basis of the knowledge of the galaxy 
distribution. We assume a King-like distribution, with param- 
eters typical of nearby/medium-redshift clusters: a core radius 
Rc = 1 /20 x 7? v i, and a slope-parameter ySfi t = 0.8, i.e. the 
volume galaxy density at large radii goes as r~ 3 ^'" = r~ 2A (see 
G98 and Girardi & Mezzetti l200"TV We obtain R P y = 1 .25 Mpc, 
with a 25% error, thus in agreement with the above direct esti- 
mate. The mass is then M(< R out ) = (1.9 + 0.5) x 1O 15 M , in 
good agreement with our first estimate. 

We can use the second of the above approaches to ob- 
tain the mass within the whole assumed virialized region, 
which is larger than that sampled by observations, M{< R v [ r = 
2.0 Mpc) = (2.2 + 0.6) x 1O 15 M . 

3. Dynamical analysis 

3.1. Velocity distribution 

We analyze the velocity distribution to look for possible devia- 
tions from Gaussianity that could provide important signatures 
of complex dynamics. For the following tests the null hypothe- 
sis is that the velocity distribution is a single Gaussian. We base 
our analysis on shape estimators, i.e. the kurtosis and the skew- 
ness. As for the kurtosis, we find K = 2.04+0.49, that indicates 
a ~ 2<r departure from a Gaussian distribution (reference value 
K = 3). In addition, we compute the scaled tail index (ST I), 
which also measures the shape of a distribution, but is based 
on order statistics of the dataset instead of its moments (see, 
e.g., Beers et al. U991> . This estimator, STI - 0.860, indicates 
that the tails are underpopulated if the parent population is re- 
ally a single Gaussian with a c.l. between ~ 90% and ~ 95%, 
(see Table 2 of Bird & Beers fT993l . Finally, also the W-test 
(Shapiro & Wilk 1965 ) rejects the null hypothesis of a Gaussian 
parent distribution at the 98% c.l.. 

Then we investigate the presence of gaps in the distribution, 
which can be the signature of subclustering. A weighted gap 
in the space of the ordered velocities is defined as the dif- 
ference between two contiguous velocities, weighted by the 
location of these velocities with respect to the middle of the 
data. We obtain values for these gaps relative to their aver- 
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Fig. 4. Velocity distribution of radial velocities for the 95 
cluster members. Bottom panel: stripe density plot where 
arrows indicate the position of significant gaps. The first gap 
lies between 249757 and 249997 km s _1 , the second be- 
tween 251886 and 252185 km s _1 , and the third between 
252305 and 252635 km s '. Top panel: velocity histogram 
with a binning of 500 km s 1 with superimposed the 
three Gaussians corresponding to KMMla, KMMlb, and 
KMM2 in Tableffl 

age size, precisely the midmean of the weighted-gap distri- 
bution. We look for normalized gaps larger than 2.25 since in 
random draws of a Gaussian distribution they arise at most in 
about 3% of the cases, independent of the sample size (Wainer 
and Schacht 1978; see also Beers et al. U991> . Three signifi- 
cant gaps (2.312, 2.366, 2.395) in the ordered velocity dataset 
are detected (see Fig.|4j. From low to high velocities the dataset 
is divided in parts containing 39, 29, 3, and 24 galaxies: thus 
the gaps individuate substantially three main subsets. 

In order to detect the presence of groups within our ve- 
locity dataset we use the Kaye's mixture model (KMM) test 
(Ashman et al. 1994 1. The KMM algorithm fits a user-specified 
number of Gaussian distributions to a dataset and assesses 
the improvement of that fit over a single Gaussian. In addi- 
tion, it provides the maximum-likelihood estimate of the un- 
known n-mode Gaussians and an assignment of objects into 
groups. KMM is most appropriate in situations where theoreti- 
cal and/or empirical arguments indicate that a Gaussian model 
is reasonable. This is valid in the case of cluster velocity dis- 
tributions, where gravitational interactions drive the system to- 
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ward a Gaussian distribution. However, one of the major un- 
certainties of this method is the optimal choice of the number 
of groups for the partition. Moreover, only in mixture mod- 
els with equal covariance matrices for all components the al- 
gorithm converges, while this is not always true for the het- 
eroscedastic case (see Ashman et al. 1994, for further details). 

Our search for significant gaps suggests the presence of two 
Gaussians (separated by the two very close second and third 
gaps at ~ 252 x 10 3 km s _1 ) or possibly three Gaussians (cor- 
responding to the three main subsets). In the homoscedastic 
case the KMM algorithm fits a two-group partition by reject- 
ing the single Gaussian at the 97.4% c.l. (as obtained from the 
likelihood ratio test). The three-group partition is fitted at the 
97.9% c.l.. In the heteroscedastic case we use the results of 
the gap analysis to determine the first guess and we fit two ve- 
locity groups around the guess mean-velocities of 249 x 10 3 
and 254 x 10 3 km s _1 . The algorithm fits a two-group parti- 
tion at the 99.4% c.l. Similarly, we fit three velocity groups 
around the guess mean- velocities of 247 x 10 3 , 250 x 10 3 , and 
254 x 10 3 km s _1 to obtain a three-group partition at the 97.2% 
c.l. The high probability value obtained in the heteroscedas- 
tic bimodal case suggests the presence of a main cluster of 76 
galaxies (KMM1), with the presence of a high-velocity clump 
of 19 galaxies (KMM2). In turn, the main cluster can be sub- 
divided in two clumps of 19 and 57 galaxies according to the 
heteroscedastic trimodal case (KMMla and KMMlb, respec- 
tively). Table^lists the kinematical properties of these clumps: 
the three corresponding Gaussians are displayed in Fig.|4| 
Fig- IU shows that the galaxies of the KMMla group mainly 
populate the South-West central region of the cluster. This 
spatial segregation suggests to investigate the velocity field in 
more detail. 

3.2. Velocity field 

The cluster velocity field may be influenced by the presence of 
internal substructures, possible cluster rotation, and the pres- 
ence of other structures on larger scales, such as nearby clus- 
ters, surrounding superclusters, and filaments. Each asymmetry 
effect could produce a velocity gradient in the cluster velocity 
field. 

To investigate the velocity field of RX J0152. 7-1357 wedi- 
vide galaxies in a low- and a high-velocity samples by using 
the median value of galaxy velocities V = 250626 km s _1 and 
check the difference between the spatial distributions of the two 
samples. High- and low-velocity galaxies appear segregated in 
the E-W direction (see Fig . |6j . The corresponding spatial dis- 
tributions are different at the 99.2% c.l. according to the two- 
dimensional Kolmogorov-Smirnov test (hereafter 2DKS-test; 
see Fasano & Franceschini 1 1 9871 as implemented by Press et 
al. lT9^2l . 

To estimate the direction of the velocity gradient we per- 
form a multiple linear regression fit to the observed velocities 
with respect to the galaxy positions in the plane of the sky (see 
also den Hartog & Katgert [T996l Girardi et al. fT996l . We find a 
position angle on the celestial sphere of PA = 97° ±16° (mea- 
sured counter-clock-wise from North), i.e. higher-velocity 
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Fig. 5. Spatial distribution on the sky of the 95 member galax- 
ies. Open symbols and crosses indicate galaxies assigned to 
KMM1 and KMM2 groups, respectively (see text). Squares and 
triangles indicate KMMla and KMMlb groups, respectively. 
The plot is centered on the cluster center defined in Sect. 2. 
Three circular regions, corresponding to regions of extended 
X-ray emission are indicated, too (see Fig. 1 by Demarco et al. 
120051 . 

galaxies lie in the East-South-East region of the cluster, in 
agreement with the visual impression of galaxy distribution 
in Fig. |6] To assess the significance of this velocity gradient 
we perform 1000 Monte Carlo simulations by randomly shuf- 
fling the galaxy velocities and for each simulation we deter- 
mine the coefficient of multiple determination (RC 2 , see e.g., 
NAG Fortran Workstation Handbook 1986 1. We define the sig- 
nificance of the velocity gradient as the fraction of times in 
which the RC 2 of the simulated data is smaller than the ob- 
served RC 2 . We find that the velocity gradient is significant at 
the 98.3% c.l.. 

We also analyze the central cluster region using 22 galax- 
ies within a radius of 0.4 Mpc. This choice allows us to include 
the position of both X-ray peaks and exclude the East region 
populated by higher-velocity galaxies only. We find a very sig- 
nificant (99.5%) position angle of PA = 59°^1, i.e. higher- 
velocity galaxies lie in the direction of the North-East X-ray 
clump. 

3.3. 3D substructure and detection of subclumps 

The existence of correlations between positions and velocities 
of cluster galaxies is a footprint of real substructures. Here we 
combine velocity and position information to compute the A- 
statistics devised by Dressier & Schectman ( 1988). This test is 
sensitive to spatially compact subsystems that have either an 
average velocity that differs from the cluster mean, or a ve- 
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Table 1. Results of Kinematical analysis 



Sample 



< V > 
km s" 1 



km s 



Whole system 



95 250530 + 135 1322+ 



KMM partitions 



KMM1 
KMM2 
KMM la 
KMMlb 



76 249730: 
19 253634: 
19 247290: 
57 250498: 



124 
50 
97 
102 



1080+ 



+ 113 

-53 

408 + ^ 



Dressier - Schectman structures 



DS - SW* 
DS -Eh- 
DS -SW 
DS -E 
DS -M 



10 

9 
76 



248338 : 
253488 : 
248613 : 
253429 : 
250510: 



146 

241 
107 
298 
149 



318!* 

545+263 

317!^ 
848!f 8 ° 
1293!% 



Properties of X - ray - centered clumps 



SW - clump(< 0.2Mpc) 
SW-clump(< 0.18Mpc) 
E - clump(< 0.4Mpc) 
NE - clump(< 0.4Mpc) 



10 

8 
7 
15 



248535 : 
248713 : 
253506 : 
251346: 



172 
121 
304 
241 



503 +4:w 

3or 122 
jui _ 107 

71 q+287 
' Iu -117 

888! 7 f 



Segregation analysis 



passive galaxies 56 250313 ±171 

"active" galaxies 39 250803 + 226 

SW - clump(passive gals)(< 0.3Mpc) 10 248488 + 1 10 



TTT 

Sl 
139 
-125 
32J + I32 
- >il -59 



1268 
1410 



o 

U 

CO 




X (arcmin) 



We use the biweigth and the gapper estimators by Beers et al. 
(1990) for samples with N g > 15 and with N g < 15 galaxies, 
respectively (see also Girardi et al. 119931 . 



"Tig. 6. Spatial distribution on the sky of the 95 member galax- 
ies: the larger the triangle, the smaller is the radial velocity. 
Open and solid triangles indicate galaxies with velocity lower 
and higher than the median cluster velocity, respectively. The 
=plot is centered on the cluster center. The big and the small ar- 
rows indicate the position angle of the cluster gradient as mea- 
sured over the whole cluster and in internal regions, respec- 
tively. The three circles correspond to the regions of extended 
X-ray emission. 



locity dispersion that differs from the global one, or both. We 
find A = 154 for the value of the parameter which gives the 
cumulative deviation of the local kinematical parameters (ve- 
locity and velocity dispersion) from the global cluster parame- 
ters. To compute the significance of substructure we run 1000 
Monte Carlo simulations, randomly shuffling the galaxy veloc- 
ities, and obtain a value of > 99.9%. 

This technique also provides information on the positions 
of substructures. Fig. Q shows the distribution on the sky of 
all galaxies, each marked by a circle: the larger the circle, the 
larger the deviation 5, of the local parameters from the global 
cluster parameters, i.e. the higher the evidence for substructure. 
A clump of galaxies with low velocity is the likely cause of 
large values of 6/ in the region which lies closely at South- 
West of the cluster center, i.e. in correspondence of the South- 
West X-ray peak. The other possible substructure, populated 
by high-velocity galaxies, lies in the Eastern region. 

To assign galaxies to the 3D-subclumps, we resort to the 
technique developed by Biviano et al. (2002 1, who used the in- 
dividual (5,-values of the Dressier & Schectman method. The 
critical point is to determine the value of 5, that optimally sep- 
arates between internal and external substructures. To this aim 
we consider the (5,-values of all 1000 Monte Carlo simulations 
already used to determine the significance of the substructure 
(see above). The resulting distribution of 5, is compared to the 
observed one finding a difference of 99.8% c.l. according to the 



KS-test. The "simulated" distribution is normalized to produce 
the observed number of galaxies and compared to the observed 
distribution in Fig. [8] the latter shows a tail at large values. 
This tail is populated by galaxies that presumably are in sub- 
structures. 

For the selection of galaxies within substructures we 
choose the value of = 3.35, since only after the rejection 
of the values > d>*, the observed and simulated distributions 
are no longer distinguishable according to the KS-test. With 
this choice, 14 galaxies of the cluster are assigned to substruc- 
tures: six to the central South- West clump (DS-SW*) and eight 
to the East clump (DS-E*), see Fig.|9] The velocity dispersions 
computed for these structures, <x v - 300 km s and =2 650 
km s _I for DS-S* and DS-E* clumps, respectively, are likely 
to be considered as lower limits since our analysis does not 
guarantee the detection of all substructure members. 

We consider also a more relaxed criteria, by selecting 
galaxies with 5j > 3 as suggested by the histogram of Fig.|8j 
Table^shows that the results for the South-West clump (DS- 
SW vs. DS-SW*) are very robust, while only one additional 
galaxy in the Eastern clump (DS-E vs. DS-E*) leads to an in- 
crease of 200 km s _1 in the velocity dispersion. We also con- 
sider the remaining 76 galaxies of the main structure (DS-M). 
DS-M does not contain significant structure according to the 
Dressler-Schectman test. However, since we cannot exclude a 
residual contamination from substructure members, its value of 
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Fig. 7. Spatial distribution of cluster members, each marked by 
a circle: the larger the circle, the larger is the deviation 6/ of the 
local parameters from the global cluster parameters, i.e. there 
is more evidence for substructure (according to the Dressier & 
Schectman test, see text). The plot is centered on the cluster 
center. 

velocity dispersion cr v ~ 1300 km s _I is likely to be an upper 
limit. 

The Dressler-Schectman results superseed those of the 
KMM test. Again, there is the presence of a low-velocity 
clump and now its South-West position is better defined by 
the detection of the DS-SW clump. The presence of a high- 
velocity clump is confirmed and located at the East by the 
detection of the DS-E clump. Moreover, the location of DS 
clumps well coincide with X-ray peaks of extended emissions. 

3.4. Analysis of X-ray — centered clumps 

The good spatial agreement between detected galaxy clumps 
and peaks of X-ray emission prompts us to analyze the pro- 
files of mean velocity and velocity dispersion of galaxy sys- 
tems corresponding to the South-West, East, and North-East 
X-ray peaks, i.e. using the position of the X-ray peaks as 
system-centers (see Figs.EB Fig. ^2 and Fig. [21 respectively). 
This allows an independent analysis of the individual galaxy 
clumps. An increase of the velocity-dispersion profile in their 
central regions might be due to dynamical friction and galaxy 
merging (e.g., Menci & Fusco-Femiano 1996; Girardi et al. 
119981 Biviano & Katgert 2004 1, or simply induced by the pres- 
ence of interlopers or of a secondary clump (e.g., Girardi et al. 
I1996> . The latter hypothesis can be investigated by looking at 
the behaviour of the mean velocity profile. Figs. [H)] and 
[T2l show velocity-dispersion and mean-velocity profiles, and 
regions likely not contaminated by other clumps and thus reli- 
able for kinematical analysis. Detailed results of this analysis 




Fig. 8. The distribution of Si deviations of the Dressler- 
Schectman analysis for the 95 member galaxies. The solid line 
represents the observations, the dashed line the distribution for 
the galaxies of simulated clusters, normalized to the observed 
number. 

are included in Table[2where the clumps are named as SW, E, 
and NE. 

The analysis of the South-West central region has indicated 
the presence of a low-velocity clump with a low velocity- 
dispersion (of 300^1-00 km s according to DS-SW and 
KMMla results). Fig. 1101 shows how the velocity-dispersion 
increases with the distance from the South-West X-ray peak. 
The mean-velocity shows a sharp change very close to the X- 
ray peak, at ~ 0.2 Mpc. This suggests a strong contamination 
of galaxies from other structures. We consider two possible 
uncontaminated regions: one within 0.2 Mpc, where we find 
cr v - 500 km s _1 , and one within 0.18 Mpc, where we find 
o" V 300kms _1 . Such a sharp change of <x v is induced just by 
the rejection of two galaxies, one of which has an anomalously 
high velocity. The value of <x v for the SW-clump is further an- 
alyzed in Sect. 3.5 and discussed in Sect. 4.1. 

The Dressler-Schectman analysis of the East region has 
indicated the presence of a high-velocity clump with a ve- 
locity dispersion of about 600-800 km s _1 . By choosing the 
X-ray peak as center (Fig. 1111 . the mean velocity changes at 
~ 0.4 - 0.5 Mpc from the X-ray peak. Inside this region, we 
obtain cr v ~ 700 km s for the E-clump 

Fig. [prefers to the region around the North-East X-ray 
peak. The main mass clump is located in this same position, 
according to the gravitational lensing analysis (Jee et al. 2005 1. 
We have shown that this region is mostly populated by galax- 
ies having velocities intermediate between those of the above 
clumps, and likely forms a high velocity dispersion struc- 
ture (i.e., KMMlb clump in Sect. 3.1, and DS-M system in 
Sect. 3.3). Fig. H2lshows an increase of the integral velocity- 




Fig. 9. Spatial distribution on the sky of the 95 member 
galaxies. Squares and crosses indicate galaxies assigned to 
the South-West and East clumps detected by the Dressler- 
Schectman analysis, respectively (DS-SW and DS-E): thick 
symbols indicate DS-SW* and DS-E*. Triangles indicate the 
remaining galaxies of the main system (DS-M). 



dispersion profile at about 0.4 Mpc from the X-ray peak, and 
a corresponding sharp change in the mean velocity. Moreover, 
galaxies of both DS-SW and DS-E substructures lie beyond 
0.4 Mpc from the North-East X-ray peak. Thus, the value 
cr v =i 900 km s _1 , computed within 0.4 Mpc , should be reli- 
able. 

The three X-ray clumps differ from each other in mean ve- 
locities at a c.l. > 99%, according to the means-test (e.g., Press 
et al. rP?92l . 

Assuming that each of the three galaxy clumps is a sys- 
tem in dynamical equilibrium, for each clump we compute the 
virial radius and the mass contained inside with the same pro- 
cedure adopted in Sect. 2 (see Table|2j. The large uncertainties 
associated to the mass values are due to poor number statistics. 

3.5. Spectral-type segregation 

We check for possible spectral-type segregation of galaxies, 
both in position and in velocity space, by using the classifi- 
cation of Demarco et al. J2005I see their Table 4), i.e. pas- 
sive galaxies (k), galaxies with significant Balmer lines - likely 
post-starbursts (k+a/a+k) and galaxies with relevant emission 
lines (e/k+a+[OII]). The sample of cluster members contains 
56, 7, and 32 passive, post-starburst, and emission-line galax- 
ies, respectively. 

Figure ^] shows the spatial distribution of galaxies of dif- 
ferent types. As already noted by Demarco et al., emission-line 
galaxies avoid the regions of the subclumps (see also Homeier 



Fig. 10. Kinematical profiles of the S W clump obtained assum- 
ing the X-ray peak as center. The vertical line indicates the re- 
gion likely not contaminated from other clumps (see Sect. 3.4). 
The dashed vertical line indicates the radius of the extended X- 
ray emission, as defined by Demarco et al. (2005). Top panel: 
rest-frame velocity vs. projected distance from the clump cen- 
ter: squares and crosses indicate the DS-SW and DS-E as in 
Fig. [5] Differential (big circles) and integral (small points) 
mean velocity and LOS velocity-dispersion profiles are shown 
in middle and bottom panels, respectively. For the differential 
profiles are plotted: a) the values for eight annuli from the cen- 
ter of the clumps, each of 0.2 Mpc (heavy symbols); the cur- 
rent values of each ten galaxies (faint symbols). For the integral 
profiles, the mean and dispersion at a given (projected) radius 
from the clump-center is estimated by considering all galax- 
ies within that radius. The error bands at the 68% c.l. are also 
shown. 



et al. 2005 1. The same behaviour is shown by post-starburst 
galaxies. When comparing spatial distributions of passive (k) 
and "active" (k+a/a+k/e/k+a+[OII]) galaxies we find a very 
strong difference: > 99.99%, according to the 2DKS-test. 

As for the velocity distributions, no difference is found be- 
tween passive and "active" galaxies, according to the KS-tests. 
Moreover, mean velocities and velocity dispersions of the two 
populations (see Tabled do not significantly differ according 
to the means- and F-test (e.g., Press et al. 1992 1. This suggests 
that our sample of member galaxies is not significantly con- 
taminated by interlopers. In fact, possible field galaxies would 
preferably contaminate the sample of "active" galaxies causing 
a difference in the kinematical properties with respect to the 
sample of passive galaxies, e.g., enhancing the velocity disper- 
sion or changing the mean velocity. 

Finally, we perform again the analysis of mean velocity and 
velocity-dispersion profiles of Sect. 3.4, by considering passive 
galaxies only. We draw different conclusions only for the SW- 
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Table 2. Virial mass estimates 
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Fig. 11. The same as in Fig. EH but referring to the E- 
clump. Note: here the annuli for the differential profiles are 
0.45 Mpc each. 
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Fig. 12. The same as in Fig. ^3 but referring to the central 
NE-clump. 



clump. Fig. 1 141 shows that the mean velocity now changes only 
at ~ 0.3 Mpc from the South- West X-ray peak and the velocity 
dispersion does not increase any longer in the central region. 
Within ~ 0.3 Mpc, we compute for the SW-clump a value of 
32 1^ 2 km s , in good agreement with the lower estimate of 
cr v obtained in Sect. 3.4 (see TableQ. 



Sample 



R yk M(<R m ) 
Mpc 10 14 xM Q 



Whole system 
SW - clump" 
SW - clump b 
E - clump 
NE - clump 



2.0 22 
0.8 
0.5 
1.1 
1.3 



1 2 +21 

-0.5 

27 +0 - 23 



3 5 +3 ° 

—1.5 

7 +2g 

'•"-2.1 



a Using the cr v computed within 0.2 Mpc (see TableQ. 
b Using the cr v computed within 0.18 Mpc (see Table0. 
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Fig. 13. Projected distribution of the 95 member galaxies. 
Circles, squares, and crosses indicate passive, post-starburst, 
and emission-lines galaxies, respectively. The three clumps 
analyzed in Sect. 3.4 are indicated by the three circles and 
correspond to the regions likely not contaminated by galaxies 
of other clumps, with radii corresponding to vertical lines in 
Figs. 10, 11, and 12. 



4. Discussion 

Out of 187 galaxies with available redshift we assign 95 mem- 
bers to RXJ0152. 7-1357. This galaxy selection is more re- 



strictive than made by Demarco et al. (2005 102 members giv- 
ing a velocity dispersion of ~ 1600 km s _1 ) due to our analysis 
of position and velocity combined together. In particular, we 
reject a small group of galaxies at z = 0.864-0.867. In spite of 
this more restrictive member selection, the value we obtain for 
the LOS velocity dispersion is still rather high, cry = 1322^ 
km s _1 , and lies in the high-tail of the cr v -distribution of 
nearby/medium redshift clusters (see, e.g., Fadda et al. 1996 
Mazure et al. 119961 Girardi & Mezzetti 12001 1 . The position 
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Fig. 14. The same as in Fig. EH but considering passive galax- 
ies only. Here the annuli for the differential profiles are 0.3 
Mpceach. 

on the Lx.t>oi-o"v plane is consistent with the relation provided 
by Borgani et al. (1999) for moderately distant clusters and 
by Wu et al. (1999) for local clusters. As for the cry-Tx re- 
lation, assuming the density-energy equipartition between gas 
and galaxies, i.e. j8 S pec = 1 (e.g., Girardi et al. 119961 ["1998 
Xue & Wu l2000b . where /3 spec = cr v /(kT/fmi p ) with pi = 0.58 
the mean molecular weight and m p the proton mass, our value 
of cry corresponds to kT = 10.6^['j keV. This value is more 
than 2<x higher than the X-ray temperature determined from 
BeppoSAX observations (Delia Ceca et al. 2000 1 and more 
than 3<x higher than those of the North-East and South- West 
X-ray systems as determined from Chandra data (Maughan et 
al. |5003, Huo et al. 2004 1. This suggests a strong departure 
from dynamical equilibrium and, in fact, we find evidence for 
non-Gaussianity of the velocity distribution, presence of a ve- 
locity gradient and significant substructure. 

We find no kinematical difference between passive and 
"active" galaxy populations. This suggests that our sample of 
member galaxies is not significantly contaminated by interlop- 
ers. In fact, possible field galaxies would preferably contami- 
nate the sample of "active" galaxies causing a difference in the 
kinematical properties with respect to the sample of passive 
galaxies, e.g., enhancing the velocity dispersion or changing 
the mean velocity. 

Instead, our analysis shows that the high value of cry is due 
to the complex structure of this system, i.e. to the presence of 
three galaxy clumps of different mean-velocity. Using optical 
data only we detect the low-velocity SW-clump in the central 
regions and the high-velocity E-clump, which lie close to the 
South- West and East peaks detected by the X-ray analysis. The 
North-East X-ray peak is then associated to the main galaxy 
structure. In particular, the high relative velocity between the 



NE- and SW-clumps, V, = (V NE - V sw )/(1 + <z» = 1531 
km s , explains the high value of cry measured in the cen- 
tral cluster region and the presence of a velocity gradient there 
(see Figs.|51and|6jl, while the global velocity gradient is induced 
by the presence of the high-velocity E-clump in external clus- 
ter regions. The presence of the three galaxy clumps was al- 
ready suggested by Demarco et al. (2005) from the inspection 
of the velocity distribution in relation to the spatial location of 
galaxies. Moreover, the NE-, SW-, and E- clumps correspond 
to three clumps in the mass distribution as obtained from the 
weak lensing analysis (Jee et al. 2005: C, F, and A subclumps, 
respectively). 

As for the mass of the whole cluster, from the global analy- 
sis of Sect. 2 we obtain M(< 2.0 Mpc) = (2.2 + 0.6) x 10 15 M o . 
Since the system is not virialized, but likely bound (see the 
discussion below), this estimate might overestimate the mass 
even by a factor two. Adding the mass estimates of each clump 
within its virial radius (see Table [2] Sect. 3.4), we obtain 
M = 1-2^Q3 x 1O 15 M : this estimate should be considered as a 
lower value within 2.0 Mpc , since it does not consider other 
small clumps or isolated infalling, bound galaxies, as well the 
likely possibility that the three clumps extend outside the virial 
radius. Thus, we conclude that the cluster mass within 2 Mpc 
lies in the range of (1.2 - 2.2) x 10 15 M Q , which is typical to 
that of very massive clusters (e.g., Girardi et al. 119981 Girardi 
& Mezzetti 20(H). 

Our mass estimate is consistent with that of Maughan et 
20031 of (1.1 +0.2)x 1O 15 M , 



al. ( 120031 of (1.1 ± 0.2) x 1O 1S M , based on Chandra X-ray 
analysis and considering only the two central clumps within 
1 .4 Mpc. To compare with results from weak lensing analyses 
we also compute the projected mass, by considering the global 
cluster geometry as formed by the three clumps at the clus- 
ter redshift. Each clump is described by the King-like mass 
distribution (see Sect. 2) or, alternatively, by a NFW profile 
where the mass-dependent concentration parameter is taken 
from Navarro et al. ©971 and rescaled by the factor 1+z 
(Bullock et al. 1200 II Dolag et al. l2004> . The mass distribution 
of each clump is truncated at one virial radius or, alternatively, 
at two virial radii. In the following we indicate the range of 
our results. We find the projected mass within 1 Mpc from the 
center of the main clump (NE-clump) to be (9 - 15)xlO 14 M , 
higher than that, 5xlO 14 M , of Jee et al. (2005), but in agree- 
ment with the value > 1 x 1O 15 M of Huo et al. (120041 see 
their Figure 10). Indeed, both Huo et al. and Jee et al. com- 
pare their weak lensing results with an isothermal sphere with 
cry =900-1000 km s , in agreement with the value of cry that 
we measure for the main galaxy clump. However, the weak- 
lensing mass lies above or below the isothermal sphere mass 
for Jee et al. and Huo et al., respectively. 



4.1. Individual clumps 

Our estimate of cry for the NE-clump well agrees with that of 
Demarco et al. (2005) and corresponds to kT = 4.8^4 keV, 
in agreement with the observed gas temperature of ~ 6 KeV 
(Maughan et al. I20031 Tozzi et al. l2003l Huo et al. 120041 
Similarly, our mass estimate, M(< i? v j r = 1.35 Mpc) = 
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7.0^2 jXl0 14 M Q , well agrees with the X-ray mass by Maughan 
et al 120031 M{< 1.4 Mpc) = 7.O+|\?xlO 14 M ]. To com- 
pare our results with other studies we rescale M(< R^) at 
their radii by using the King-like profile or, alternatively, the 
NFW profile (see above). In the following we give the two val- 
ues obtained from the rescaling, reliable with a 30% lower- 
error and a 40% upper-error, as derived from the estimate 
of M(< R vu ). Our estimates well agree with those of other 
studies: M(< 0.43 Mpc) =(1.9-2.6)xlO 14 M , cf. with (2.5 + 
O.9)xlO 14 M by Demarco et al. (2005 1, based on galaxy dy- 
namics; M(r < 65 "= 0.496 Mpc) =(2.3-2.9)xl0 14 M G , cf. with 
(3 + l)xlO 14 M by Joy et al. ( |2"0"0"T1 . based on the Sunyaev- 
Zeldovich effect; M(r < 0.753 Mpc)=(3.8- 4.3)xl O 14 M , cf. 
with (2.66 + O.77)xlO 14 M by Ettori et al. dlUOll . based on 
Chandra X-ray data; M(r < 1 Mpc) =(5.1-5.4)xlO 14 M , cf. 
with ~ 5xlO 14 M by Huo et al. (2004i, based on Chandra X- 
ray data. 

As for the SW-clump, the results in the literature are not yet 
clear. In fact, the X-ray temperature suggests that the North- 
East and the South-West clumps are similar in mass (Maughan 
et al. 120031 Huo et al. I2004> . while both the optical analysis 
by Demarco et al. ( 2005 1 and the weak lensing analysis by Jee 
et al. ( 2005 1 find that the South-West clump is about half mas- 
sive than the North-East clump. Our analysis of the cr v -profile 
gives two alternative values for cry: the larger value is consis- 
tent with that found by Demarco et al. ( 2005 cf. cry = 503+ 4 g 9 
km s _1 vs. their 737 + 126 km s~') and with the observed 
gas temperatures of 5-6 keV (Maughan et al. 2003; Huo et al. 
2004i, while the lower estimate, cry = 301^[ ( p, is significantly 
different. This uncertainty is due to the fact that the <x v pro- 
file increases in central regions (see FigllOl and thus the cry 
estimate strongly depends on the considered region. Demarco 
et al. considered a region (based on X-ray data) larger than 
our region (based on kinematical data). Our analysis of pas- 
sive galaxies also gives a small value, cry ~ 300 km s _1 , 
thus suggesting two alternative hypothesis: 1) high values of 
cry are due to galaxy-contamination by other clumps, so that 
the SW-clump should be considered as a very small group, 
2) we are dealing with a very relaxed core hosted in a high- 
cry, massive cluster. The second hypothesis is consistent with 
the observations of nearby clusters where cry of the subsam- 
ple of bright central elliptical galaxies is lower than cry of the 
whole cluster (Biviano & Katgert 2004), a phenomenon possi- 
bly due to dynamical friction and galaxy merging (e.g., Menci 
& Fusco-Femiano 1996 1. Only a deeper galaxy sample would 
allow us to better trace and separate the North-East and the 
South-West systems and thus discriminate between the two hy- 
potheses. However, the SW-clump appears to be so dense of 
galaxies that we are inclined to believe in the detection of a 
cluster-core. In this case, we note that: a) our mass estimate 
would be an underestimate of the global mass of the Southern 
cluster; b) our results would be reconciled with high values of 
gas temperature and X-ray luminosity (Maughan et al. 2003 
Tozzi et al. l2U0H Huo et al. l2U04l ). 

As for the Eastern clump, the level of X-ray emission in 
the Chandra image is much lower than those of the North-East 
or the South-West clumps (see Fig. 1 of Demarco et al. 2005 1. 
On the contrary, its gravitational-lensing mass is comparable 



to that of the South-West clump (see A and F clumps by Jee 
et al. 2005 1, and our estimate of velocity dispersion is typical 
of a massive cluster, cry ~ 700 km s . This discrepancy with 
X-ray luminosity suggests that this galaxy system is far from 
being virialized, maybe elongated along the LOS (thus giving a 
high cr v and a high projected lensing mass), with the gas com- 
ponent not very dense. In particular, the Eastern X-ray peak 
might be associated to a small group embedded in a large- 
scale structure filament connecting to the cluster from the 
East-South-East region, which is populated by higher ve- 
locity -maybe more distant - galaxies (see Fig. In the 
case of a bound, but non virialized structure, we might have 
overestimated the mass even by a factor two. 

We finally compare the projected mass of the three clumps 
within a radius of 20 " with the results from weak lensing by Jee 
et al. (2005 1. The resulting values for projected masses of the 
NE-, SW-, and E-clumps lie in the ranges (1.6 - 2.3)xlO 14 M , 
(0.5 - O.7)xlO 14 M , and (1.0 - 1.5)xlO 14 M , all values some- 
what higher than those reported in Table2 of Jee et al. for 
clumps C, F, and A, respectively. 

4.2. Analytic calculations of the dynamical status 

Here, we investigate the relative dynamics of the NE- and SW- 
clumps in the central cluster region using different analytic ap- 
proaches which are based on an energy integral formalism in 
the framework of locally flat spacetime and Newtonian grav- 
ity (e.g., Beers et al. 119821 . The values of the relevant observ- 
able quantities for the two-clump system are: the relative LOS 
velocity, V r = 1531 km s , the projected linear distance be- 
tween the two clumps, D = 0.66 Mpc, the mass of the system 
obtained adding the masses of the two clumps each within its 
virial radius, M sys = 8.2+^xl0 14 M G . 

First, we consider the Newtonian criterion for gravita- 
tional binding stated in terms of the observables as V 2 D < 
2GM sys sin 2 acosa, where a is the projection angle between the 
plane of the sky and the line connecting the centers of two 
clumps. The faint curve in Fig.^]separates the bound and un- 
bound regions according to the Newtonian criterion (above and 
below the curve, respectively). Considering the value of M sys , 
the NE+SW system is bound between 30° and 77°: the cor- 
responding probability, computed considering the solid angles 

rll 

(i.e., J 3 „ cos a da), is 47%. We also consider the implemented 
criterion V 2 D < 2GMsin 2 aycosao, which introduces differ- 
ent angles for projection of distance and velocity, not assuming 
strictly radial motion between the clumps (Hughes et al. 1995 1. 
We obtain a binding probability of 44%. 

Then, we apply the analytical two-body model introduced 
by Beers et al. i 19821 1 and Thompson (119821 see also Lubin 
et al. 1998 for a recent application). This model assumes ra- 
dial orbits for the clumps, with no shear or net rotation of the 
system. Furthermore, the clumps are assumed to start their evo- 
lution at time to = with separation do = 0, and are moving 
apart or coming together for the first time in their history, i.e. 
we are assuming that we are seeing the cluster prior to merging. 
The bimodal model solution gives the total system mass M sys 
as a function of a (e.g., Gregory & Thompson 1984). Fig. 1151 
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Fig. 15. System mass vs. projection angle for bound and un- 
bound solutions of the two-body model applied to the NE- and 
SW-clumps (solid and dotted curves, respectively, see text). 
The thin curve separates the bound and unbound regions ac- 
cording to the Newtonian criterion (above and below the curve, 
respectively). The horizontal lines give the observational values 
of the mass system and its lcr error bands. 

compares the bimodal-model solutions with the observed mass 
of the system, which is the most uncertain observational pa- 
rameter. The present bound outgoing solutions (i.e. expanding), 
BO, are clearly inconsistent with the observed mass. The pos- 
sible solutions span these cases: the bound and present incom- 
ing solution (i.e. collapsing), BI a and Bit, and the unbound- 
outgoing solution, U O. For the incoming case there are two so- 
lutions because of the ambiguity in the projection angle a. We 
compute the probabilities associated to each solution assum- 
ing that the region of M sys values between lcr bands is equally 
probable for individual solutions: f BIa ~ 65%, P B ib ~ 35%, 
Puo < 1 x 10~ 4 %. 

There are several limitations to characterize the dynamics 
of the central region of RXJ0152. 7-1357 using these mod- 
els. For instance, possible underestimates of the masses, e.g., 
if the clumps extend outside the virial radius or if the SW- 
clump is only the core of the South-West system (see above), 
lead to binding probabilities larger than those computed above. 
Moreover, the models do not take into accounts the mass dis- 
tribution in the clumps when the separation of the clumps is 
comparable with their size (i.e. at small a) and do not con- 
sider the possible effect of the E-clump. Finally, the two-body 
model breaks down in a regime where the NE- and SW-clumps 
are already strongly interacting, as suggested by several evi- 
dences: the displacement between peaks of gas distribution and 
of galaxy/dark matter distribution (Maughan et al. 2003 , Huo 
et al. l2004> : Jee et al.l2"005): the possible presence of a shock 
front (Maughan et al. 2003 1; the presence of galaxies showing 



a very recent star formation episode (J0rgensen et al. 2005 1; 
the segregation of star-forming and non star-forming galax- 
ies probably induced by the interaction with the intra-cluster 
medium (Homeier et al. l2005> . 

Looking at galaxies only we cannot discriminate between a 
pre- or post-merging phase since the galaxy component is very 
robust against mergers, e.g., two clusters can pass through one 
another without destroying the individual optical components 
(e.g., White & Fabian [T9951 Roettiger et al. lT997l . Note, for 
instance, that the properties of the SW-clump resemble those 
of cluster-cores destined to survive tidal disruption during the 
merger: size comparable to the cluster core and mass < 0.05 x 
cluster mass (see Gonzalez-Casado et al. 1994 1. These cores 
will be detectable in the host cluster as a substructure for a long 
time. Since the gas component shows two well distinct entities 
in the central cluster regions, we presume that the merging is 
not too advanced, i.e. well before the coalescence. 

Under the assumption that the two central clumps are al- 
ready very close, we apply the above dynamical models to the 
system made of the [(NE+SW)+E] clumps, too. The values 
of the relevant observable quantities are: V, = 1401 km s _1 , 
D = 1.09 Mpc, and M sys = H-^'l xK) I4 M . We obtain 
that the binding probabilities are 48%, and 45%, according to 
the Newtonian criterion and its implementation, respectively; 
while the two-body model gives a probability > 99.9% for the 
bound incoming solution. 

5. Summary & conclusions 

We present the results of the dynamical analysis of the cluster 
of galaxies RXJ0152. 7-1357, one of the most massive struc- 
tures known at z > 0.8. The X-ray emission is known to have 
two clumps in the central regions, and a third clump ~ 1 Mpc 
to the East. Our analysis is based on velocities and positions of 
member galaxies taken from the extensive spectroscopic survey 
performed by Demarco et al. (120051 . i.e. 187 galaxies having 
redshift in the cluster region. 

We find that RXJ0152. 7-1357 appears as a well isolated 
peak in the redshift space at z = 0.836, and select 95 cluster 
members. We compute a value for the LOS velocity dispersion 
of galaxies, <x v = 1322^ km s _1 , much larger than expected 
for a relaxed cluster with an observed X-ray temperature of 
-5-6 keV. 

We find evidence that this cluster is far from dynamical 
equilibrium, as shown by: 

- the non Gaussianity of the velocity distribution according 
to different tests, at the 90-98% c.l., as well as the presence 
of significant velocity gaps; 

- the correlation between velocities and positions of galaxies 
at the > 99% c.l., and the presence of a velocity gradient; 

- the presence of significant substructures at the > 99.9% c.l.. 

To detect and analyze possible subsystems we used differ- 
ent methods. 

- By applying the KMM method we find that a two-clumps, 
and likely a three-clumps partition of the velocity distri- 
bution is significantly better than a single Gaussian to de- 
scribe the velocity distribution; in particular, the galaxies of 
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KMMla group are mainly located in the South-West cen- 
tral region. 

- By combining positions and velocities in the Dressier & 
Schectman statistics we detect two substructures, well cor- 
responding in location to the South-West and East X-ray 
peaks, in addition to the main cluster component identified 
with the North-East X-ray peak. 

- Taking advantage of X-ray peak determination, we analyze 
the three galaxy clumps centered in these peaks through the 
profiles of mean velocity and velocity dispersion. This anal- 
ysis allows us to estimate the clump region that is likely not 
contaminated by galaxies of other clumps and to evaluate 
the kinematical properties. 

In summary, our analysis shows that the high value of 
<t v is due to the complex structure of RXJ0152. 7-1357, i.e. 
to the presence of three galaxy clumps of different mean- 
velocity. Using optical data we detect a low-velocity clump 
(with cr v =300-500 km s -1 ) in the central South-West re- 
gion and a high-velocity clump (with cr v ^700 km s -1 ) in 
the Eastern region, nicely matching the position of the South- 
West and East peaks detected in the X-ray emission. The cen- 
tral North-East X-ray peak is associated to the main galaxy 
structure having intermediate velocity and cr v ~ 900 km s . 
The three clumps differ from each other in mean velocities at a 
c.l. > 99% (relative LOS velocities are > 1000 km s -1 ). 

The mass of the whole system within 2 Mpc is estimated to 
be (1.2-2.2)xl0 15 M o , where the upper and lower limits come 
from the virial analysis of the cluster as a whole and from the 
sum of virial masses of the three individual clumps, respec- 
tively. 

Analytic calculations, based on the two-body model, indi- 
cate that the system is most likely bound, destined to merge. In 
particular, we suggest that the South-West clump is not a small 
group, but rather the dense core of a massive cluster, able to 
survive tidal disruption during the merger. 

In conclusion, RXJ0152. 7-1357 reveals a very complex 
structure, with several clumps likely destined to merge in a 
very massive cluster. Our results lend further support to the pic- 
ture that massive clusters at z > 0.8 are dynamically complex 
and, therefore, likely to be young. This indicates that we are 
approaching the epoch at which such massive structures take 
shapes from the evolution of the cosmic web. On-going exten- 
sive spectroscopic surveys of such systems at z ~ 1 and beyond, 
combined with detailed analyses of their gaseus and dark mat- 
ter components (now possible with weak lensing analysis of 
HST-ACS data; Jee et al. l2005l Lombardi et al. l2005l . will shed 
new light on cluster formation processes. 
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